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ABSTRACT 

Galaxy and Mass Assembly (GAMA) is a project to study galaxy formation and evolution, 
combining imaging data from ultraviolet to radio with spectroscopic data from the AAOmega 
spectrograph on the Anglo-Australian Telescope. Using data from phase 1 of GAMA, taken 
over three observing seasons, and correcting for various minor sources of incompleteness, we 
calculate galaxy luminosity functions (LFs) and their evolution in the ugriz passbands. 

At low redshift, z < 0.1, we find that blue galaxies, defined according to a magnitude- 
dependent but non-evolving colour cut, are reasonably well fitted over a range of more than 
ten magnitudes by simple Schechter functions in all bands. Red galaxies, and the combined 
blue-plus-red sample, require double power-law Schechter functions to fit a dip in their LF 
faintwards of the characteristic magnitude M* before a steepening faint end. This upturn is at 
least partly due to dust-reddened disc galaxies. 

We measure evolution of the galaxy LF over the redshift range 0.002 < z < 0.5 both by 
using a parametric fit and by measuring binned LFs in redshift slices. The characteristic lumi- 
nosity L* is found to increase with redshift in all bands, with red galaxies showing stronger 
luminosity evolution than blue galaxies. The comoving number density of blue galaxies in- 
creases with redshift, while that of red galaxies decreases, consistent with prevailing move- 
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ment from blue cloud to red sequence. As well as being more numerous at higher redshift, 
blue galaxies also dominate the overall luminosity density beyond redshifts z ~ 0.2. At lower 
redshifts, the luminosity density is dominated by red galaxies in the riz bands, by blue galax- 
ies in u and g. 

Key words: galaxies: evolution — galaxies: luminosity function, mass function — galaxies: 
statistics. 



1 INTRODUCTION 

Measurements of the galaxy luminosity function (LF) and its evolu- 
tion provide important constraints on theories of galaxy formation 
and evolution, (see e.g. |Benson et al.|2003| l. It is currently believed 
that galaxies formed hierarchically from the merger of sub-clumps. 
Looking back in time with increasing redshift, the star formation 
rate appears to peak at redshift z ~ 1, above which it plateaus and 
slowly declines towards z ~ 6 (|Cole et al.||2000| |Hopkins|20fJ4) 



Hopki ns & Beacom|2006l |Yuksel et d.|2008[|Kistier et al.|2009f 
Since z — 1, galaxies are thought to have evolved mostly pas- 
sively as their stellar populations age, with occasional activity trig- 
gered by accretion and interactions with other galaxies. |Noeske et] 
[aL] ( |2007 i) have suggested that the first major burst of star formation 
is delayed to later times for low mass galaxies, contributing to the 
downsizing phenomenon. 

There has long been a discrepancy between the measured 
number density of low-luminosity galaxies (the 'faint end' of the 
LF) and predictions from cold dark matter (CDM) hierarchical sim- 
ulations, in the sense that fewer low-luminosity galaxies than pre- 
dicted by most models are observed (Trentham & Tully 20021. Of 
course, interpretation of these simulation results is subject to un- 
certainties in the baryon physics. In particular, more effective feed- 
back in low mass halos might act to suppress the faint end of the 
LF. However, it is also possible that many surveys have underesti- 
mated the number of dwarf galaxies due to the correlation between 
luminosity and surface brightness which makes them hard to de 



tect (Driver 1999] |Cross & Driver|2002l|Cameron & Driver|2007l 



|2009| >. |G"eiler et al.| ( [20TT) have recently demonstrated that the LF 
faint-end slope steepens with decreasing surface brightness. 

Galaxy LFs have previously been measured in the ugriz 
bands from the Sloan Digital Sky Survey ( SDSS,|York et al.|2000> 
by Blanton et al. 1 2003b)> , [Loveday |{2004|>,|Blanton et al.H2" 005b), 
Montero-Dorta & Prad"a|j2009| >, and |Hill et al.|p010| >. |Blanton et| 
aL] 1200 3b) analysed a sample of 147,986 galaxies, roughly equiv- 



alent to SDSS Data Release 1 (DRl, [Abazajian et al.|2003) . They 
fit the LF with a series of overlapping Gaussian functions, allowing 
the amplitude of each Gaussian to vary, along with two parameters 
Q and P describing, respectively, luminosity and density evolu- 
tion. They maximized the joint likelihood of absolute magnitude 
and redshift, rather than the likelihood of absolute magnitude given 
redshift, making this estimator more sensitive to evolution, as well 
as to density fluctuations due to large-scale structure. They found 
luminosity densities at z = 0.1 to increase systematically with ef- 
fective wavelength of survey band, and for luminosity evolution 
to decline systematically with wavelength. Allowing for LF evo- 
lution enabled reconciliation of previously discrepant luminosity 
densities obtained from SDSS commissioning data (Blanton et al. 
200 1|| an d the Two-degree field Galaxy Redshift Survey l Folkes et 



al.|1999||Norberg et al.|2002) . 



galaxies brighter than M r — 5lgh = —21.5 mag was a factor ~ 3 
higher at redshift z = 0.3 than today, due to luminosity and/or 
density evolution. 

|Blanton et a l. ( 2005b I focused on the faint end of the LF of 
low-redshift galaxies from SDSS DR2 ( Abazajian et al. 2004}, and 
found that a double-power-law Schechter function was required to 
fit an upturn in the LF at M r — 51g/i> — 18 mag with faint- 
end slope «2 — —1.5 after correcting for low surface-brightness 
incompleteness. 

|Montero-Dorta & Prada| ([2009) have analysed SDSS DR6 
l |Adelman-McCarthy et al.|2008^ which is ro ughly five times larger 
than the sample analysed by Bl anton et al.| ( |2003b) . Their results 
are generally consistent with those of Blanton et al., although they 
do point out a bright-end excess above Schechter function fits, par- 
ticularly in the u and g bands, due primarily to active galactic nuclei 
(AGNs). A bright-end excess above a best-fitting Schechter func- 
tion has also been observed in near-IR passbands by lones et al. 
((2006). 



Hill et al.] ( |2010p analysed combined datasets from the Millen- 

— , SDSS and the UKIDSS 



nium Galaxy Catalogue i Liske et al 



2003 



Large Area Survey ( |Lawrence et al. 



2007 



over a common volume 
0. 1 to obtain LFs in the 



of ~ 71, 000/i~ 3 Mpc a within redshift z 
ugrizY JHK bands. They found that LFs in all bands were rea- 
sonably well fitted by Schechter functions, apart from tentative up- 
turns at the faint ends of the i- and z-band LFs. Hill et al. provided 
the first homogeneous measurement of the luminosity density (LD) 
over the optical-near-IR regimes, finding a smooth spectral energy 
distribution (SED). 

Here we present an estimate of ugriz galaxy LFs from the 
Galaxy and Mass Assembly (GAMA, |Driver et al.1[2009| [20TT) 
survey. GAMA provides an ideal sample with which to constrain 
the galaxy LF at low to moderate redshifts due to its combination 
of moderately deep spectroscopic magnitude limit (r < 19.4 or 
r < 19.8) and wide-area sky coverage (three 4 x 12 deg 2 regions). 

We describe the input galaxy sample and incompleteness, ve- 
locity and A'-corrections in Section [2] Our LF estimation proce- 
dure is described in Section [3] and tested using simulations in Ap- 
pendix]^ We present our results and a discussion of luminosity and 
density evolution in Section]?] with our conclusions summarized in 
Section|5] 

Unless otherwise stated, we assume a Hubble constant of 
Ho = 100ft km/s/Mpc and an Qm = 0.3, Qa = 0.7 cosmology in 
calculating distances, co-moving volumes and luminosities. 



2 DATA AND OBSERVATIONS 
2.1 Input catalogue 



[Loveday (2004) measured the r-band LF in redshift slices 
from SDSS DR1 and found that the comoving number density of 



The input catalogue for GAMA is described in detail by B aldry et] 
[ah (2010 1. In brief, it consists of three 4 x 12 deg 2 regions cen- 
tred approximately on the equator and at right ascensions of 9, 12 
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and 14.5 hours. These fields are known as G09, G12 and G15 re- 
spectively. Primary galaxy targets were selected from Data Release 
6 (DR6, |Adelman -McCarth y et al. 12008} of the Sloan Digital Sky 
Survey (SPSS, | York et al.]20 00) to extinction-corrected, Petrosian 
magnitude limits of r < 19.4 mag in the G09 and G15 fields and 
r < 19.8 mag in the G12 field. 

We require Petrosian and model magnitudes and their errors in 
all five SDSS passbands in order to determine Jf-corrections (Sec- 
tion [23}, and so we match objects in the GAMA team catalogue 
TilingCatvl6 to objects in the SDSS DR6 PhotoObj table on SDSS 
ObjID using the SDSS CASJOB^] service. We use only objects 
with GAMA SURVEY_CLASS > 3 in order to exclude additional 
filler targets from the sample. We exclude objects, which, upon vi- 
sual inspection, showed no evidence of galaxy light, were not the 
main part of a deblended galaxy, or had compromised photometry 
(VIS-CLASS = 2, 3 or 4 respectively). See |Baldry et al.| l [20To} for 
further details of these target flags and Section |2.7| for a discussion 
of additional visual inspection of extreme luminosity objects. 

In estimating LFs, we use Petrosian magnitudes corrected 
for Galactic extinction according to the dust maps of Schlegel, 
|Finkbeiner, & Da vis ( 19981. We make no attempt here to correct 
for intrinsic dust extinction within each galaxy, as was done by 
|Driver et al. ( 2007}, nor to extrapolate the Petrosian magnitudes to 
total, as done, for example, by |Graham et al.||2005[ > and [Hill et al.| 
< |201 1| >. These systematic corrections to SDSS photometry, much 
more significant than any small random errors, will be considered 
in a subsequent paper. 

An exception to our use of Petrosian magnitudes is for u-band 
data, where we instead use a pseudo-Petrosian magnitude defined 
by 

^pseudo — Pctro — ^modcl fmodcl ~t~ ^"pctro- (1) 

The reason for this is that the Petrosian u-band quantities are noisy 
and suffer from systematic sky-subtraction errors (Bald ry et al.| 
[2005}. The pseudo-Petrosian u-band magnitude defined above, (us- 
ing the SDSS r band since it has highest signal-to-noise), and re- 
ferred to AS lisclcct by |Baldry et al.|p005| >, is much better behaved 
at faint magnitudes. 

For colour selection (see Section |2~6} , we use SDSS model 
magnitudes in defining (g — r) colour, as recommended by the 
SDSS websiteE] 

2.2 Spectroscopic observations 

GAMA spectroscopic observations are described in the first 
GAMA data release paper (Driver et al. 2011). Observations for 
the GAMA Phase 1 campaign were made over 100 nights between 
2008 February and 2010 May, comprising 493 overlapping two- 
degree fields. Redshifts were assigned in a semi-automated fashion 
by the observers at the telescope. A subsequent re-redshifting ex- 
ercise (Liske et al. in prep.) was used to assign a normalised qual- 
ity nQ to each redshift, according to each particular observer and 
their assigned quality Q. Here we use reliable (nQ > 2) redshifts 
from all three years of the GAMA Phase 1 campaign. In addi- 
tion to pre-existing redshifts and those obtained with the Anglo- 
Australian Telescope, twenty redshifts of brighter galaxies were 
obtained with the Liverpool Telescope. The GAMA-II campaign, 
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extending the survey to additional southern fields, began in 2011, 
but only GAMA-I redshifts are used here. 



2.3 Completeness 

Although GAMA has a very high spectroscopic completeness (> 
98 per cent; Driver et al. 201 1), the small level of incompleteness is 
likely to preferentially affect low surface brightness, low luminos- 
ity galaxies, or galaxies lacking distinctive spectral features. We 
have identified three sources of incompleteness that potentially af- 
fect the survey: the completeness of the input catalogue (imaging 
completeness), completeness of the targets for which spectra have 
been obtained (target completeness) and the success rate of ob- 
taining spectroscopic redshifts (spectroscopic success rate). These 
three sources of incompleteness, and how we correct for them, are 
now considered in turn. 



2.3.1 Imaging completeness 

Imaging completeness has been estimated for the SDSS main 
galaxy sample by Bla nton et al.H2005b} , who passed fake galaxy 
images through the SDSS photometric pipeline. Blanton et al. 
found that imaging completeness is nearly independent of appar- 
ent magnitude (at least down to m r ~ 18 mag), depending mostly 
on r-band half-light surface brightness, /xso, r (their Fig. 2). Thus, 
while GAMA goes about 2 mag fainter than the SDSS main galaxy 
sample, the Blanton et al. imaging completeness should still be ap- 
proximately applicable. We have used their imaging completeness 
estimates modified in the following waysQ 

(i) Blanton et al. determine imaging completeness over the sur- 
face brightness range 18 < /iso,r < 24.5 mag arcsec -2 . Extrap- 
olating their completeness as faint as fi 50 ^ r — 26 mag arcsec -2 
results in negative completeness values. We therefore arbitrarily as- 
sume 1 per cent imaging completeness at /iso.r = 26 mag arcsec -2 
and linearly interpolate from the faintest tabulated Blanton et al. 
completeness point at (/iso,r, / P h) = (24.34, 0.33). 

(ii) The Blanton et al. imaging completeness decreases at the 
bright end, /iso.r ^ 19 mag arcsec -2 , due to a lower angular size 
limit of 6*50 > 2 arcsec and a star-galaxy separation criterion 
A sg = r ps { — r mo dd > 0.24 for the SDSS main galaxy sam- 
ple which excludes some compact, high surface-brightness (HSB) 
galaxies. GAMA target selection uses a far less stringent A sg > 
0.05, backed up by J — K colour selection, and so is much more 
complete in HSB galaxies. We therefore omit the Blanton et al. 
completeness points at a 5(l>r < 19.2 mag arcsec -2 and instead as- 
sume 100 per cent completeness at ^so.r = 19.0 mag arcsec -2 and 
brighter. 

Our revised imaging completeness curve, along with a histogram 
of Mso.r values for GAMA galaxies, is given in Fig[T| 



3 An alternative way of estimating imaging completeness is to determine 
what fraction of galaxies detected in the much deeper co-added data from 
SDSS Stripe 82 ( Abazaji anet al.|2009> are detected in regular, single-epoch 
SDSS imaging. However, one needs to allow for the large number of bright 
star or noise images misclassified as low surface-brightness galaxies in the 
SDSS co-added catalogue, and so this approach was abandoned. It will be 
re-explored once high-quality VST KIDS imaging of the GAMA regions 
becomes available. 
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Figure 2. GAMA target completeness as a function of magnitude (pseudo-Petrosian for u, Petrosian for griz). The upper panels show galaxy counts in 
varying-width magnitude bins, chosen to give roughly equal numbers of galaxies per bin, for GAMA targets (thin black histogram) and counts for galaxies 
that have been spectroscopically observed (thick red histogram). The lower panels show target completeness, ie. the ratio of observed to target counts, with the 
horizontal dotted line indicating 100 per cent target completeness. Vertical dashed lines indicate our chosen magnitude limits in each band: 20.3, 19.5, 19.4, 
18.7, 18.2 in ugriz respectively. 




Pa, r / mag arcsec 2 



Figure 1. Top line shows imaging completeness as a function of r-band 
half-light surface brightness, li5o,r, from Blanton et af| {200 5b), modi- 
fied as described in the text. The histogram shows the normalised counts 
of ^50, r for galaxies in the GAMA sample. 



2.3.2 Target completeness 

Target completeness in the r band may be assessed relative to the 
GAMA tiling catalogue, which contains all galaxies to r — 19.8 
mag in the GAMA regions. In the ugiz bands, however, there is 
no well-defined magnitude limit. We therefore re-implement the 
GAMA target selection criteria detailed by Baldry et al. (2 010| > on 
samples of objects selected from the SDSS DR6 PHOTOOBJ table. 
We replace the Baldry et alT] l |2010) magnitude limits (their equa- 



tion 6) with the following: u < 21.0, g < 20.5, r < 20.0, i < 19.5 
or z < 19.0. 

Target completeness in each band is then simply defined as 
the fraction of target galaxies that have been spectroscopically ob- 
served, either by GAMA or by another redshift survey, as a func- 
tion of apparent magnitude in that band. This is shown in Fig. [2] 
where we have used magnitude bins which are equally spaced in 
m! — I0 0,52 ( m-mmin '. This binning is chosen to give a roughly 
equal number of galaxies per bin, thus avoiding large Poisson un- 
certainties at bright magnitudes. In the r band, target completeness 
is around 98-99 per cent brighter than r = 19.4 mag corresponding 
to the magnitude limit of the GAMA G09 and G15 fields. 

In the other four bands, the drop in completeness at faint mag- 
nitudes is more gradual due to the spread in galaxy colours. Mag- 
nitude limits in each band are set to the faintest magnitude at which 
target completeness is at least 92 per cent (it band) or where com- 
pleteness starts to drop rapidly. These magnitude limits are 20.3, 
19.5, 19.4, 18.7, 18.2 in ugriz respectively and are indicated by 
the vertical dashed lines in Fig. [2] 

An alternative approach to estimating the LF in bands other 
than that of target selection is to perform a multivariate LF, e.g. 
|Loveday| ( |2000| >, or to use a 1/V max estimator where V ma j( is cal- 
culated using the selection-band magnitude, e.g. Smit h, Loveday &] 
[Crossl|2009l . While using more data, these estimators suffer from 
a colour bias as one approaches the sample flux limit, and so we 
prefer the more conservative approach adopted here. 

2.3.3 Redshift success rate 

Redshift success rate is most likely to depend on the flux that goes 
down a 2dF fibre, that is a seeing-convolved 2-arcsec-diameter 
aperture. The closest quantity available in the SDSS database is 
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Figure 3. GAMA redshift success rate as a function of fibre r-band magni- 
tude plotted as a black histogram, with the horizontal dotted line indicating 
100 per cent success. The red curve shows the best-fit sigmoid function. 

f iberMag_r, hereinafter rahic, corresponding to the flux con- 
tained within a 3-arcsec-diameter aperture centred on the galaxy. 
We therefore determine histograms of rgbrc (uncorrected for Galac- 
tic extinction) for all objects with high-quality redshifts (nQ > 2) 
and for all objects with spectra. The ratio of these two histograms 
then gives redshift success as a function of rfibrc, and is shown in 
Fig. [3] Note that some spectra observed in poor conditions have 
been re-observed at a later date in order to obtain this high success 
rate. 

We see that redshift success rate is essentially 100 per cent for 
?"fibrc < 19.5, declines gently to around 98 per cent by rfibrc = 20 
and then declines steeply at fainter magnitudes. We have fitted a 
sigmoid function / = 1/(1 + e a ^ rribTe ~ b ') to the binned success 
rate. Sigmoid functions have previously been used to model survey 
completeness, e.g. by[Ellis & Bland-Hawthorn ( 2007 ). Our best-fit 
sigmoid function has parameters a = 1.89 mag - , b = 21.91 mag, 
shown by the red line in Fig. [5] and we use this fit in determining 
redshift success rate. 

2.3.4 Galaxy weights 

Each galaxy is given a weight which is equal to the reciprocal of the 
product of the imaging completeness, target completeness and red- 
shift success rate. Imaging completeness Ci m is determined from 
the galaxy's apparent r-band half-light surface brightness, ^50, r by 
linear interpolation of the curve in Fig[T] Target completeness Ctarg 
is determined separately in each band from the galaxy's magnitude 
according to Fig. [2] and the spectroscopic success rate C spec is de- 
termined from the sigmoid function fit described in section |2.3.3| 
The weight Wi assigned to galaxy i is then 

Wi = l/(Ci m CtargC s pec)- (2) 

These weights, as a function of magnitude in each band, are shown 
for a randomly selected 5 per cent of galaxies in Fig. [4] The ma- 
jority of galaxies brighter than our magnitude limits have weight 
Wi < 1.1, with a small fraction extending to Wi — 2. 

2.4 Velocity corrections 

The redshifting software RUNZ (Saunders, Can non & Sutherland! 
|2004| > provides heliocentric redshifts. Before converting heliocen- 
tric redshifts to any other velocity reference frame, we first elimi- 



nate likely stars from our sample by rejecting objects with a helio- 
centric redshift Zhciio < 0.002 (cz < 600 km/s). This lower red- 
shift cut is conservatively chosen, as the 2nd Catalogue of Radial 
Velocities with Astrometric Data ( Kharchenko et al. 2007 1 includes 
only one star with radial velocity RV > 500 km/s, and the vast ma- 
jority of Galactic stars have RV < 200 km/s. Furthermore, Fig. 7 
of |Driver et al.| ( [201 1[ > shows that the redshift-error distribution for 
GAMA is essentially zero by cz = 600 km/s. 

Having eliminated 2111 likely stars from our sample, he- 
liocentric redshifts are converted to the CMB rest frame zcmb 
according to the dipole of Linewea ver et al.| (T996). For nearby 
galaxies (.zcmb < 0.03), we apply the multiattractor flow model 
of |Tonry et al.| ( |2000) . Note there are triple-valued solutions of 
zcmb — > ZTonry over a small range in G12 (near the Virgo cluster), 
here the average distance is used. The solution is tapered to zcmb 
from zcmb = 0.02 to zcmb = 0.03 (see |Baldry et al.|20"TT| for de- 
tails). We will see later that the Tonry et al. flow correction affects 
only the very faintest end of the LF. 



2.5 A'-corrections 

When estimating intrinsic galaxy luminosities, it is necessary to 
correct for the fact that a fixed observed passband corresponds 
to a different range of wavelengths in the rest frames of galaxies 
at different redshifts. the so-called A-correction (Humason, May- 
all & Sandage 1956). The if -correction depends on the passband 
used, the redshift of the galaxy and its spectral energy distribu- 
tion (SED). Here we use KCORRECT v4_2 (Blanton et al.|2003al 
|Blanton &"R oweis 2007) in order to estimate and apply these cor- 
rections. Briefly, this code estimates an SED for each galaxy by 
finding the non-negative, linear combination of five template spec- 
tra that gives the best fit to the five SDSS model magnitudes of 
that galaxy. KCORRECT can then calculate the if -correction in any 
given passband at any redshift. Before calling KCORRECT itself, we 
use K_SDSSFIX to convert SDSS asinh magnitudes to the AB sys- 
tem and to add in quadrature to the random ugriz mag errors given 
in the SDSS database typical systematic errors of (0.05, 0.02, 0.02, 
0.02, 0.03) mag respectively. 

We determine if -corrections in a passband blueshifted by 
zo = 0.1. These magnitudes are indicated with a superscript prefix 
of 0.1, e.g. 01 M r . This choice of rest frame allows direct compar- 
ison with most previously published LFs based on SDSS data. 

Our LF estimators require A'-corrections to be calculated for 
each galaxy at many different redshifts in order to work out visibil- 
ity limits. To make this calculation more efficient, we fit a fourth- 
order polynomial Pk (z) = Ysi=o a i( z ~ z o) % , withzn = 0.1, to the 
K(z) distribution for each galaxy and use this polynomial fit to de- 
termine A'-corrections as a function of redshift. Using a polynomial 
of this order, the rms difference between the KCORRECT prediction 
and the polynomial fit is 0.01 mag or less in all five bands. Calcu- 
lated A'-corrections and their polynomial fits are shown for the first 
four galaxies in our sample, along with the median A'-corrections 
and the 5 and 95 percentile ranges for the full sample, in Fig. [5] 

Strictly, one should use heliocentric redshift when calculat- 
ing A'-corrections, since they depend on the observed passband. 
However, for consistency with finding the minimum and maximum 
redshifts at which a galaxy would still be visible when using the 
1/Knax LF estimator, we use the flow-corrected redshift as de- 
scribed in section |2~4] The difference in A'-correction when using 
heliocentric or flow-corrected redshift is entirely negligible. 
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Figure 4. Completeness-correction weights as a function of magnitude for a random 5 per cent subset of GAMA galaxies. The vertical dashed lines show the 
magnitude limits applied in the LF analysis of each band. 



2.6 Colour sub-samples 

As well as analysing flux-limited samples of galaxies in the ugriz 
bands (hereafter the combined sample), we separate the galaxies 
into blue and red sub-samples. Following Zehavi et al. (20fTJ, we 
use a colour cut based on 7f -corrected (g — r) model colour and 
absolute r-band magnitude that is insensitive to redshift: 

0,1 (g - r) mod ei = 0.15 - 0.03 1 M r , (3) 

We have adjusted the zero-point of 0.21 mag in |Zehavi et al.lpoTT) 
to 0.15 mag in order to better follow the 'green valley' and to get 
more equal-sized samples of blue and red galaxies. This colour cut 
works well at all redshifts (Fig.[6|, although we see that the colour 
bimodality becomes far less obvious beyond redshift z — 0.2 due 
to the lack of low-luminosity, blue galaxies at these high redshifts. 

Although colour bimodality is more pronounced in (u — r) 
colour, e.g. Strat eva et akl pOOl I, Bald ry et al1 < |2004} , u-band pho- 
tometry, even after forming a 'pseudo-Petrosian' magnitude (equa- 
tion [TJ is rather noisy, and so we prefer to base our colour cuts 
on the more robust (g — r) colour. This colour cut (in the original 
form of Zehavi et al.) has also been used to investigate the angular 
clustering of galaxies by Christodoulou et aL] < |201 1 p . One should 
also note that colour is not a proxy for galaxy morphology: many 
red galaxies are in fact dust-obscured disc galaxies (Driver et al. in 
prep.). 

2.7 Outlier inspection 

We measure the LF over a very wide range of luminosities, —23 to 
— 1 1 mag in the r band. Galaxies at the extremes of this luminos- 
ity range are very rare in a flux-limited survey, due either to their 
intrinsic low number density at high luminosity or small detection 
volume at low luminosity, and thus it is likely that a significant 
fraction of these putative extreme objects are in fact artifacts due to 
incorrectly assigned redshifts or magnitudes. The first author has 



Table 1. Classification of extreme high- and low-luminosity objects. 



POSTXLASS 01 M u <-20mag 01 M r >-15mag 

1 OK 4,743 299 

2QSO 18 

3 Major shred 68 62 

4 Minor shred 7 

5 Problem deblend 151 16 

6 Bad sky background 246 14 



therefore inspected image cutouts showing SDSS image detections 
of 5,226 very luminous GAMA targets with 0,1 M u < -20 mag 
and 398 very faint targets with °'°M r > —15 mag. We choose the 
u band to select luminous outliers since the u-band LF shows the 
largest bright-end excess. 

Table[T]shows how the inspected images were classified. The 
classification codes, which we call POST-CLASS in order to distin- 
guish them from the pre-target-selection VIS_CLASS have the fol- 
lowing meanings. 

(1) OK — nothing from the image would lead one to expect 
poor photometry for that object. 

(2) The object looks like a QSO, i.e. blue and point-like. This 
classification is ignored in the analysis (treated as OK) due to the 
difficulty of distinguishing QSOs and compact blue galaxies from 
the imaging data alone. 

(3) The central part of a galaxy which has been shredded into 
multiple detections. It is likely that the luminosity is somewhat un- 
derestimated in these cases. 

(4) The target is a minor part of a shredded galaxy. Luminosity 
is likely to be greatly underestimated. 

(5) The galaxy is very close to a second object which either has 
not been deblended, or is likely to significantly affect the estimated 
luminosity in either direction. 
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Figure 5. Top four panels: /^-corrections as a function of redshift (red, blue, 
green, magenta and yellow respectively for ugriz) for the first four galaxies 
in our sample. Black dashed lines show fourth-order polynomial fits to each 
band. Bottom panel: median ^-corrections (coloured continuous lines) and 
5 and 95 percentile ranges (dotted lines) for the entire GAMA sample. Black 
dashed lines show fourth-order polynomial fits to the medians. 
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Figure 6. 0,1 (g — r) colour versus 01 M r r-band absolute magnitude con- 
tour plots for GAMA galaxies in four redshift ranges as labelled. Ten con- 
tours, spaced linearly in density, are colour-coded from black to red in or- 
der of increasing density. The straight line shows the magnitude-dependent 
colour cut separating blue and red galaxies given by equation (3). 



(6) Photometry is likely severely compromised by rapidly vary- 
ing sky background due typically to the presence of a nearby satu- 
rated star. 

Examples of objects with these classifications are shown in Fig. [7] 
In practice, there is some ambiguity in assigning a galaxy with clas- 
sification 4 or 5, but as far as the LF analysis is concerned, it makes 
no difference. 

These inspections were based on version 10 of the GAMA 
tiling catalogue, excluding objects with VIS_CLASS = 2—4. The ma- 
jor and minor shreds in Table[T]have also been inspected by 1KB. In 
the case of major shreds, we have summed the flux from the compo- 
nents of the shredded galaxy to derive a 'deblend-fixed' magnitude. 
In total, 281 GAMA-I galaxies have had their magnitudes fixed in 
this manner. 

In the case of six minor shreds which both JL and 1KB agreed 
on, the value of VIS_CLASS has been set equal to 3 in the latest 
version (vl6) of the tiling catalogue. 

The fractions of galaxies with POST_CLASS of 3 or 4 or higher 
as a function of M u and M r are shown in Fig. [8] We see in the left 
panel that by magnitudes of M u = —20, less than 10 per cent of 
objects have suspect photometry. Once we allow for fixing of over- 
deblended galaxies, a similar fraction of objects with M r ~ — 15 
have suspect photometry (right panel). 

For our analysis, we have chosen to exclude any galaxies with 
POST_CLASS of 4 or higher, i.e. we include major shreds with fixed 
fluxes but exclude minor shreds, problem deblends and bad sky 
objects. Fig. [9] shows the ratio of the r and u band LFs using 
POST-CLASS < 4 galaxies to that determined using all galaxies. 
We see that excluding objects with suspect photometry has a rela- 
tively minor effect on the determined LF: the very bright end and 
some faint-end bins are systematically lower by up to 50 per cent; 
these changes are comparable to the size of the error bars. 

Finally, we note that Brough et al. (201 1 1 have independently 
checked a sample of GAMA galaxies with the lowest detected Ha 
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Figure 7. Examples of objects classified from 1 to 6. The GAMA targets are 
at the centre of each image, which are 40 arcsec on each side. The colour 
table has been inverted, so that red objects appear blue and vice-versa, in 
order to obtain a light background. Circles denote SDSS image detections: 
multiple circles on a single object (example classifications 3 and 4) suggest 
that it has been over-deblended. POST.CLASS classifications are shown in 
the top-left corner of each image, their meaning is given in Table[T] 
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Figure 8. Fraction of objects with POST_CLASS > 2 (blue circles) and > 3 
(green squares) as a function of 01 M u (left panel) and 0,1 M r (right panel). 
Error bars show Poisson errors on the counts. 
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Figure 9. Ratio of LFs determined using POST.CLASS < 3 objects to that 
using all objects in r band (top) and u band (bottom). Symbols show ratio 
of SWML estimates and their uncertainties, the continuous lines shows the 
ratio of parametric fits to the two samples. 



luminosity. Our four faintest r-band luminosity galaxies are also in 
the B rough et al. sample. 



3 ESTIMATING THE LUMINOSITY FUNCTION AND 
ITS EVOLUTION 

3.1 Parameterizing the evolution 

In order to parametrize the evolution of the galaxy LF, we follow 
|Lin et al.|(l999) in assuming a SchechtendT976b function in which 
the characteristic magnitude M* and galaxy density 0* are allowed 
to vary with redshift, but where the faint-end slope a is assumed to 
be non-evolving FJ 

Specifically, in magnitudes, the Schechter function is given by 



4>{M) =0.41nlO?i*(10 ' 4(M 
where the Schechter parameters a, M* and 



"))!+«« e X p(-10°' 4(M *- M) ), 

(4) 

vary with redshift 



a(z) = a(z ), 
M*{z) =M*(z )-Q{z-zo), 



(5) 



Here the fiducial redshift zq is the same redshift to which magni- 
tudes are if-corrected, in our case Zq = 0.1. The Schechter param- 
eters a, M*(zo) and (f>*(0) and evolution parameters Q and P are 
determined via the maximum-likelihood methods described by Lin 
|etaLl ( fl999) . 

4 Evolution in the LF faint-end slope a with redshift is still rather poorly 
constrained. Ellis et al. 1 1996 1 claim that a steepens with redshift, due to an 
increase in the number of faint, star-forming galaxies at z ~ 0.5. jllbert et| 



alj(|2005 


also measure a 


Liu et al. 


{2008} find that 



(2008) find that a gets shallower at higher redshifts. Our assump- 



tion of fixed a is largely based on practical necessity, since a can only be 
reliably measured at redshifts z < 0.2 from the GAMA data. 
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First, the shape parameters a, M*{zq) and luminosity evolu- 
tion parameter Q are fit simultaneously and independently of the 
other parameters by maximising the log likelihood 



ln£ = W.lnp,. 



(6) 



Here, Wi is the incompleteness correction weighting (equation |2j 
and the probability of galaxy i having absolute magnitude Mi given 
its redshift z t is 

//rnin[M max (2^ ) ,M 2 ] 
/ <t>(M)dM , 

(7) 

where Mi, M2 are the absolute magnitude limits of the sample, 
Af m i n (zi), M max (zi) are the minimum and maximum absolute 
magnitudes visible at redshift Zi, and (f>{M) is the differential LF 
given by {4|. 

The density evolution parameter P and normalization <j>*(0) 
cancel in the ratio in ([7J and so must be determined separately. Lin 
et al. show that the parameter P may be determined by maximising 
the second likelihood 



(8) 



where 



p' t =p[zi\Mi(0),Q] 



= 10 u 



min[z max [M i (0),z 2 ] 



W u ^ z ^-dz, (9) 
az 



where z\, Z2 are the redshift limits of the sample, Zmi n , z max are 
the redshift limits over which galaxy i may be observed, given 
the survey's apparent magnitude limits and its absolute magnitude 
evolution-corrected to redshift zero, Mi(0) = Mi(zi) + Qzi. Note 
that the value of P is independent of the fiducial redshift zq. 

Finally, we fit for the overall normalisation 4>*(0). We depart 
slightly from the prescription of |Lin et aT] fT999 ) here. In place of 
their equation 14, we use a minimum variance estimate of the space 
density n of galaxies: 



E 



WiUjzj) 



dV 

dz 



S(z)U(z)dz , (10) 



where S(z) is the galaxy selection function, U(z) a redshift 
weighting function chosen to give minimum variance and dV/dz is 
the volume element at redshift z. The selection function for galax- 
ies with luminosities L\ to L2 is 

/"min(L max (2),L 2 ) / rL 2 

S{z)= / <j>(L,z)dL I I (j>(L,z)dL. (11) 

J max(L m j n (2) ,Li ) / J L\ 

Note that the integration limits in the numerator depend on the as- 
sumed AT-correction. In this case, we use the median /^-correction 
of the galaxies in the sample under consideration: see Fig[5]for me- 
dian /^-corrections for the full sample as a function of redshift. Our 
results change by much less than the estimated 1-sigma errors (see 
Section[33J if we use mean instead of median /^-corrections. 
We adopt the redshift weighting function 



U(z) = 



l + 4ir(n/W)J 3 (r c )S(z)' 



Mrc) = 



2 £{r)dr, 
(12) 



where £(r) is the two point galaxy correlation function and W is 



the mean incompleteness weighting. Provided Ja(r c ) converges 
on a scale r c much smaller than the depth of the survey, then 
this redshift weighting sch eme (equation |1 2} mini mizes the vari- 
ance in the estimate of n (Davis & Huchra 19821. Larger values 
of J3 weight galaxies at high redshift more highly; we assume 
4-7T J3 = 30, 000/i -3 Mpc 3 . This value comes from integrating the 
flux-lim ited two-point galaxy correlation function of Zeha vi et al.| 



(2005 



CM = (r/5.59)" 1S4 ,tor c = 60/i _1 Mpc; at larger sepa- 
rations the value of J3 becomes uncertain. However, the results are 
not too sensitive to the value of J3, the estimated densities changing 
by less than 8 per cent if J3 is halved. This possible systematic er- 
ror is generally comparable to or less than the statistical uncertainty 
in <f>* (5- 25 per cent). 

We check our minimum variance normalisation by comparing, 
in Tables[3]|4]and[5] the observed number of galaxies in each sam- 
ple (within the specified apparent magnitude, absolute magnitude 
and redshift limits) with the prediction 

^ /-Zmax /--tmaxfz) 



prcd 



W 



4>(L, z)-^-dz. 



(13) 



3.2 Luminosity density 

Given our assumed evolutionary model, the predicted LD is given 
by 

p Lflt = Pi (0)10°- 4(P+Q)z , (14) 
( |Lin et al.|l 999 1 equation 11), where 

p L {0) = J L<P(L, z = 0)dL = <j>*(0)L*(0)r{a + 2), (15) 

and T(x) is the standard Gamma function. In making this predic- 
tion, we are integrating over all possible luminosities, and hence 
extrapolating our Schechter function fits. This extrapolation intro- 
duces no more than 1 per cent in additional LD beyond that con- 
tained within our luminosity limits. We obtain luminosities in solar 
units using the following absolute magnitudes for the Sun in SDSS 
band passes: °''M ^ - 5 lg h = 6.80, 5.45, 4.76, 4.58, 4.51 
mag pianton et al.|2003bY 

We also directly determine LD as a function of redshift by 
summing the weighted luminosities of galaxies in a series of red- 
shift bins: 



PL 



=-y 



■ Vj ^ S L (z,)' 



(16) 



(Lin et al. equation 16). Here Vj is the volume of redshift bin j, the 
sum is over each galaxy i in bin j and the factor 

/rnin(.L max (2),L2) / roc 

Sl(z)= / L(j>(L,z)dL / / L<j}{L,z)dL, 

imax(L mill (z),li) / JO 

(17) 

(Lin et al. equation 17) extrapolates for the luminosity of galaxies 
lying outside the accessible survey flux limits. 



3.3 Binned LF estimates 

In order to assess how well the model (equation [3J parametrizes 
LF evolution, we also make non-parametric, binned, estimates of 

dSchmidt 



the LF in independent redshift ranges using the 
1968[|Eales|1993| > and stepwise maximum likelihood (SWML,[Ef 



stathiou, Ellis & Peterson 19881 methods. We use 60 magnitude 
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Figure 10. Illustration of estimating <fi(M, z) in bins of absolute magnitude 
M and redshift z represented by dotted lines for a fictitious survey with 
apparent magnitude limits m-bright = 14.5 and m,f a ; nt = 18. Galaxies 
(represented by points) are of course only found between these flux limits, 
corresponding to the lower and upper curved lines respectively. Now con- 
sider the highlighted bin, centred on M = —19.5 and with redshift limits 
z\ a = 0.10, Zfo = 0.15. At the lower redshift limit, the absolute magni- 
tude corresponding to mf alnt is M ~ —19.6. Since this is mid-bin, the 
LF estimated for this bin would be underestimated, and therefore the bin 
should be excluded. Thus for the redshift slice 0.10 < z < 0.15, only 
magnitude bins brightward of M = — 20 should be used. (The fact that the 
magnitude bin centred on M = —20.5 is incomplete at redshifts z > z\ a 
will be compensated for by 1/Vmax weighting.) A similar incompleteness 
may arise for bins at low redshift and high luminosity. For the redshift slice 
0.00 < z < 0.05, only magnitude bins fainter than M = —21 should be 
used. 



bins from M = -25 to M = -10 with AM = 0.25 and a series 
of redshift slices. 

When estimating the LF over restricted redshift ranges, one 
has to be careful to only include magnitude bins that are fully sam- 
pled, since otherwise the LF will be underestimated in incompletely 
sampled bins, see Fig. [TO] We therefore set the following magni- 
tude limits for each redshift slice so that only complete bins are 
included: 



A/faint < m&int ~ DM(z lo ) — K{z\o), 
Mbright > Tlbright — DM(z h i) — K(Zhi) 



(18) 



Here, mfaint and mbright are the flux limits of the survey, DM(z) 
is the distance modulus, K(z) is the if -correction for a galaxy with 
the median SED of those in the survey, z\ and are the limits of 
the redshift slice under consideration, and Mf a i nt and Mbright are 
the absolute magnitude limits of each bin. A bin should only be 
included if it satisfies both equations ( [18) . 

Again following |Lin et aT7| ^1999( 1, we incorporate the galaxy 
incompleteness weights into the SWML maximum likelihood esti- 
mator by multiplying each galaxy's log-probability by its weight 
before summing to form a log-likelihood (equation |6j. In the 
1/Vmax estimate, we form a sum of the weight of each galaxy di- 
vided by the volume within which it is observable. We normalize 
the SWML estimates <^>swml in each redshift slice to the 1/Vmax 
estimates cj>v max by imposing the constraint 



0SWML fc V(M ft ) = 

fc=l fe=l 



(19) 



xk V(M k ), 

where V(Mk) is the volume (within the redshift limits of each 



slice) over which a galaxy of absolute magnitude M k , being the 
mean galaxy absolute magnitude in bin k, is visible. The predicted 
number of galaxies 



w J2 c/) k V(M k )AM 



(20) 



may also be compared with the observed number of galaxies within 
each redshift range. 

We can use our SWML LF estimates to assess the quality of 
the parametric fits using a likelihood ratio test (Efstat hiou, Ellis &| 
|Peterson|1988) . In this test, we compare the log-likelihoods ln£i 
and In C2 given by equation |6j for the SWML and parametric es- 
timates respectively. The log likelihood ratio —2 Ym^Ci/Cz) is ex- 
pected to follow a x 2 distribution with v — N p — 4 degrees of 
freedom. Here N p is the number of bins in the stepwise estimate 
and we subtract 1 degree of freedom for each of the fitted shape 
parameters a, M* (0) and Q and for the arbitrary normalisation. 

To allow for the finite bin sizes and redshift ranges of the 
SWML estimates, we calculate binned estimates of the parametric 
fits. These are given by ( jLin et al.|1999[ > 



rM k +AM/2 r min[z„ 
J M k -AM/2 J max 2,, 



1 



M k +AM/2 (.min[z max (M),z 2 ] 
M k -AM/2 Jmax[z min (M), Z1 ] 



<p(M,z)^dz dM 



(21) 

Here, the parametric LF <f>(M, z) is weighted by the number 
of galaxies at each magnitude and redshift, given by the factor 



4>{M, z) ' 



These binned versions of the parametric fits are also 



used when plotting the LFs. For absolute magnitudes in all plots, 
we use the weighted mean magnitude of the galaxies in each bin, 
rather than using the magnitude of the bin centre. This helps to 
overcome the bias due to the finite width of magnitude bins. 



3.4 LF faint end 

It is now widely recognised that a Schechter function provides a 
poor fit to galaxy LFs when measured over a wide range of magni- 
tudes (e.g. Blanton et al. 2005b). In order to parametrize the faint 
end, we separately analyse a low redshift (z < 0.1) subset of the 
data and fit (non-evolving) double power-law Schechter functions. 
We use the parameterization of Loveday ( 1997 1, namely 



<f>{L) = <f>* 



exp 



—L 

17 



1 + 



(22) 



In this formulation, the standard Schechter function is multiplied by 
the factor [1 + (L/Lt)^], where L t < L* is a transition luminosity 
between two power-laws of slope a(L <C L t ) and a+/3 (L 3> Lt). 
It is fitted to unbinned data using an extension to the method of 
|Sandage, Tammann & Yahll] ( |1979^ . With this four-parameter fit 
(the normalisation <f>* is fitted separately), one has to be careful to 
choose sensible starting values in order for the downhill simplex 
algorithm (scipy . optimize . fmin) not to get stuck in local 
minima of — ln£. (We also found that it helped to call the mini- 
mizer several times, using 'best-fitting' parameters from one func- 
tion call as starting parameters for the next; — In C was found to 
converge with 2-4 calls of the minimizer.) 

Note that the double power-law Schechter function may equiv- 
alently be written as the sum of two Schechter functions, e.g. |Blan^| 
|ton et al.|p005b| >; |Baldry, Glazebrook, Driver| ( |2008) , a fact which 
comes in useful when integrating the LF. 
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Table 2. Change in fitted Schechter parameters for combined samples when 
applying imaging completeness correction. 



Band 


Aa 


AM* /mag 


A^/h'^Mpc' 3 


u 


-0.05 


-0.03 


0.00017 


9 


-0.05 


-0.04 


-0.00031 


r 


-0.06 


-0.07 


-0.00062 


i 


-0.05 


-0.06 


-0.00051 


z 


-0.03 


-0.03 


-0.00011 



When fitting a double power-law Schechter function, the like- 



lihood ratio test has 
tion|33l. 



5 degrees of freedom (cf. sec- 



3.5 Error estimates 

Schechter and evolution parameter estimates are strongly corre- 
lated, and so in Section|4]we present 95 per cent likelihood contour 
plots of shape parameters a, M* , /3, Mt and evolution parame- 
ters Q, P. For uncertainties in tabulated measurements, we esti- 
mate errors using the jackknife technique, as follows. We divide 
the GAMA survey area into nine regions, each 4x4 deg 2 . We 
then calculate the LF and LD using the methods discussed above, 
omitting each region in turn. For any parameter x, we may then 
determine its variance using 



Var(:r) = 



N 



N 



1 N 



(23) 



where N = 9 is the number of jackknife regions and x is the mean 
of the parameters Xi measured while excluding region i. The jack- 
knife method has the advantage of providing error estimates which 
include both uncertainties in the fitted parameters as well as sample 
variance. 

The sample variance in galaxy density n may also be deter- 
mined using: 



SnY 
n J 



4nJ 3 
V ' 



(24) 



Pavis & Huchra| 1982| [Efst athiou, El lis & Peterson| 1988) , where 
J3 is defined in l |12) and V is the volume of each sample between 
redshift limits. 

For errors on binned LFs, we use Poisson errors for 1/V ma x 
estimates and an inversion of the information matrix for SWML 
estimates (Efstafhiou, Ellis & Peterson 19881. 




Figure 11. LF estimates in the r band for low-redshift galaxies (2 < 0.1) 
with (solid symbols and line) and without (open symbols, dashed line) ap- 
plying a correction for imaging completeness. Symbols show SWML esti- 
mates, lines show best-fitting Schechter functions. 




-18 -16 
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Figure 12. LF estimates in the r band for low-redshift galaxies (z < 0.1) 
using the CMB reference frame (open circles) and the Tonry et al. (2000) 
multiattractor flow model (filled circles) using the SWML estimator. Solid 
and dashed lines show the best-fit Schechter functions which are indistin- 
guishable. 



4 RESULTS 

Before presenting our main results, we first check the effects of 
correcting for imaging completeness and the choice of flow model 
in converting redshifts to distances. 

4.1 Imaging completeness correction 

In Fig. [TT| we compare r-band LFs calculated for the combined 
sample, with distances calculated using the |Tonry et al.| ( (2000) mul- 
tiattractor flow model, with and without the correction for imag- 
ing completeness described in Section [2.3. 1| As expected, we see 
that applying imaging completeness corrections boosts the LF faint 
end, while barely changing the bright end. The changes in fitted 
Schechter parameters due to imaging completeness corrections are 



tabulated in Table [2] Future plots and tabulated parameters will in- 
clude imaging completeness corrections; approximate uncorrected 
Schechter parameters may be obtained by subtracting the appropri- 
ate quantities listed in Table|2] 

4.2 Effects of velocity flow model 

Luminosities of galaxies at the extreme faint end of the LF, being 
very close by, will be sensitive to peculiar velocities. In Fig. [12] 
we compare r-band LFs calculated using the CMB reference frame 
l |Lineweaver et al.|1996) and the Tonry et al. ( 2000 1 multiattractor 
flow model for galaxies at low redshift (z < 0.1). We see that the 
CMB -frame and flow model LFs only begin to differ at the extreme 
faint end, M r — 5 log ft > — 15 mag, and even at these faint mag- 
nitudes the differences are not large given the size of the error bars. 
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Figure 13. ugriz LFs at low redshift (z < 0.1). Black squares show SWML estimates for combined red and blue samples, blue circles and red triangles 
show SWML LFs for the blue and red samples respectively. Open symbols of the same shapes show the corresponding 1 /V max estimates — these are hidden 
beneath the SWML estimates for all but the very faintest galaxies. Continuous lines show the best-fit non-evolving double power-law Schechter function fits, 
dotted lines show standard Schechter function fits. LFs for the blue and red samples have been scaled by a factor of 0. 1 to aid legibility. Open diamonds show 
the 'corrected' LF (without colour selection) from Blanton et al. ( 2 005b} . 
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In particular, the recovered Schechter parameters are indistinguish- 
able between the two velocity frames. Subsequent analysis will use 
the Tonry et al. flow model to determine luminosities. 



4.3 Luminosity function faint end 

Having looked at the effects of incompleteness and flow correc- 
tions, we now study in detail the faint end of the LF for low red- 
shift (z < 0.1) galaxies. Fig.[T3]shows the LFs for our three (com- 
bined, blue and red) samples in the ugriz bands. Also shown are 
LFs corrected for surface brightness incompleteness by Bl anton et| 
[aL| (|2005b) from the New York University Value-Added Galaxy 
Catalog (NYU-VAGC) low-redshift sample < |Blanton et al.|2005a| >. 
Since these Blanton et al. LFs were calculated using restframe K- 
corrections, we apply an offset of 2.51g(l + zn) to their absolute 
magnitudes in order to convert to our zo = 0.1 band-shifted K- 
corrections. Our faint-end LFs are systematically lower than those 
of Blanton et al., particularly in the u band. The difference can 
largely be explained by the different flow models used by Blanton 
et al. and in the present analysis. Re-analysing the Blanton et al. 
data using the Tonry et al. flow model results in much better agree- 
ment JBaldry etal.|2011} — the extra 1.7 mag depth of the GAMA 
versus the SDSS main galaxy sample means that uncertainties due 
to the flow model affect the measured LF only at a correspondingly 
fainter magnitude. 

Table [3] shows the number of galaxies and absolute magni- 
tude limits for each sample, along with the parameters of stan- 
dard Schechter function fits and luminosity densities. Only for blue 
galaxies in the u, i and z bands does a standard Schechter function 
provide a statistically acceptable fit to the data at the 2 per cent level 
or better. For red galaxies, we observe a decline in number density 
faintwards of the characteristic magnitude A/* with a subsequent 
increase in faint-end slope at Mt ~ M* + 3. For the red galaxies, 
and the combined sample, a double-power-law Schechter function 
l |22fr is required to fit the shape of the observed LFs. These findings 
are in apparent agreement with the predictions of halo occupation 
distribution models, e.g. |Brown et al.| ( |2008[ >, in which luminous 
red galaxies are central galaxies, but fainter red galaxies are in- 
creasingly more likely to be satellites in relatively massive halos. 
An alternative perspective is provided by Peng et al. (2010J, who 
explain the change in faint-end slope of red galaxies via a simple 
picture for the quenching of star formation by the distinct processes 
of 'mass quenching' and 'environment quenching' . Our results pro- 
vide the most precise demonstration of the changing faint-end slope 
of red galaxies to date. 

However, the observed upturn needs to be interpreted with 
caution, since, from a quick visual inspection, the 164 faint 
( 01 M r — 51g/i > —16 mag), red galaxies that comprise the up- 
turn include a significant fraction (~ 50 per cent) of galaxies that 
appear to be disc like, as well as a number of artifacts. It thus seems 
likely that dust-reddened disc systems, as well as dwarf galaxies 
with intrinsically red stellar populations, contribute to the faint-end 
upturn in the red galaxy LF. Future work will investigate the LF 
dependence on morphology and dust reddening, utilising GAMA's 
multi-wavelength coverage. 

Double-power-law Schechter function fits are given in Ta- 
ble [4] Likelihood ratio tests show that the double-power-law 
Schechter function provides significantly better fits than the stan- 
dard Schechter function for the combined and red galaxy samples, 
at least for the redder bands. For the blue galaxies, however, the 
double-power-law Schechter function fits are actually worse than 
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Figure 14. 2-<r likelihood contours for various parameter pairs in double 
power-law Schechter function fits to the combined sample for ugriz bands 
as labelled. 



the standard Schechter function fits when taking into account the 
two additional degrees of freedom. 

The quoted errors need to be treated with caution due to strong 
correlations between the parameters, particularly in the case of the 
five-parameter, double power-law Schechter function fits. Fig. |14| 
shows 2a likelihood contours for each pair of parameters from 
a,M*,/3 and M t for the combined sample. We see that a and 
P individually are very poorly constrained, with an uncertainty 
Act ~ A/3 ~ 0.5. However, the overall faint-end slope a + /3 is 
very well constrained, with a consistent value in all five passbands 
a + /3 = — 1.37±0. 05 for the combined sample. For blue galaxies, 
Q + /3 = - 1. 50 ±0.03 and for red galaxies, a + /3 = -1.6 ±0.3. 
Consistent faint-end slopes are found for the stellar mass function 
(Baldry et al. 2011) . Also from Fig. [14] we see that the charac- 
teristic magnitude M* is positively correlated with slope a but 
negatively correlated with ft. The transition magnitude M t is only 
weakly correlated with either slope parameter a or j3 and almost 
completely uncorrelated with the characteristic magnitude M*. 

Fig. [T5] shows 2cr likelihood contours for the red sample. We 
see that, while still uncertain, the slope parameters a and j3 are only 
weakly correlated. The characteristic magnitude M* is positively 
correlated with a but almost completely uncorrelated with /3. The 
transition magnitude Mt is strongly anti-correlated with a and vir- 
tually independent of (3. 

Since the shape of the blue galaxy LF is reasonably well fitted 
by a standard Schechter function, there are huge degeneracies be- 
tween the double power-law Schechter function parameters, and so 
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Table 3. Standard Schechter function fits for low-redshift galaxies. Samples are as given in the first column. Mi and 0,1 M2 are the absolute magnitude 
limits, iV ga i the number of galaxies in the sample and N pIc ^ the predicted number of galaxies from integrating the LF (equation |l3| . a, 01 M* and <j>* are 
the usual Schechter parameters, and Pg t gives the probability of the Schechter function describing the observed LF determined from the likelihood ratio test 
described in the text. Luminosity densities pLfa ar >d PL sum are calculated via equations |15| and \16) respectively. 
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Table 4. Double power-law Schechter function fits for low-redshift galaxies. Values for 0,1 M\, 0-1 M% and -/V ga i are the same for each sample as in Table|3] 
Other columns have the same meaning as the previous Table and in addition we tabulate the values of the double power-law Schechter parameters f} and Mt- 
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the contour plots contain no useful information, and hence are not 
shown. 

In summary, in analysing the faint end of the LFs, we have 
found that: 

(i) While a standard Schechter function provides an acceptable 
fit to the blue galaxy LF in all bands, the red galaxy LF exhibits a 
decline just faintwards of M* followed by a pronounced upturn at 
magnitude Mt — M* + 3. Such an LF is well-fitted by a double- 
power-law Schechter function. 

(ii) We caution that the faint end of the red galaxy LF is possibly 
dominated by dust-reddened systems, rather than by galaxies with 
intrinsically red stellar populations. 

(iii) Neither standard nor double-power-law Schechter function 
faint-end slopes show any systematic dependence on passband: 



while strongly colour-dependent, faint-end slopes are largely in- 
dependent of passband. 

(iv) The characteristic magnitude M* (and to a lesser extent, 
the transition magnitude M t ) brightens systematically and signifi- 
cantly with passband effective wavelength. 



4.4 Luminosity function evolution 

We present LFs for the combined, blue and red samples in the 
ugriz bands in four redshift ranges in Fig. [16] Table [3] gives the 
magnitude limits (chosen to exclude the upturn seen in the LF of 
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given in Table [5] and shown below in Fig.[T9] we sub-divide into 
eight redshift bins, perform least-squares fits to the SWML esti- 
mates in each, and then fit straight lines to lg <f>* and M* versus 
redshift. We now discuss LF evolution separately, for the u, g and 
riz bands. 



4.4.1 u-band evolution 

We observe a gradual brightening of M* for all samples between 
z ~ 0.05 and z ~ 0.15, and already a bright-end excess above a 
Schechter function is becoming apparent. By z ~ 0.25, the evolv- 
ing model provides a very poor fit to the non-parametric estimates: 
the former is much steeper than the latter, and by z ~ 0.4 the 
parametric fit over-predicts the number density of galaxies by al- 
most an order of magnitude. From the least-squares fits to the non- 
parametric estimates, we see a dramatic brightening of M* in the 
highest redshift range (there are too few red galaxies to obtain a 
sensible LF fit in this bin). This is due to the very shallow slope at 
the bright end of the LF, leading the Schechter fit to prefer bright- 
ening M* to increasing <f>* . In reality, we suspect that this shallow 
slope is caused by an increasing fraction of highly-luminous AGNs, 
rather than by such strong luminosity evolution in non-active galax- 
ies. An alternative explanation is that photometric errors in the u 
band are manifesting themselves as unrealistically strong luminos- 
ity evolution. This possibility will explored when VST KIDS data 
become available in the GAMA regions. 



Figure 15. 2-cr likelihood contours for various parameter pairs in double 
power-law Schechter function fits to the red sample for ugriz bands as 
labelled. 

red galaxiesQ observed and predicted numbers of galaxies, and 
Schechter and evolution parameters in each band. Qualitatively, 
the riz LFs appear to be well fitted by the parametric evolution 
model, although this model is formally excluded by the likelihood 
ratio test for almost every colour, band and redshift combination. 
Even by eye, we see that the evolving Schechter function fits are 
in extremely poor agreement with the u and g band non-parametric 
(SWML and 1 /Vmn-x) estimates in the highest redshift range, in the 
sense that the model overpredicts the number density of luminous 
galaxies by almost an order of magnitude in the u band. Equation|5] 
thus provides a poor fit to evolution of the it- and g-band LFs be- 
yond z ~ 0.2 and z ~ 0.3, respectively. It is very possible that 
the M-band flux of more luminous, higher redshift galaxies is be- 
ing dominated by AGN. We intend to investigate the LFs of AGN- 
dominated/starforming/quiescent galaxies in a future paper. 

In addition to the parametric fit, we also fit Schechter functions 
to the SWML estimates for each redshift slice using least-squares. 
Because the LF faint end is poorly sampled at redshift z > 0.1, we 
only fit for all three Schechter parameters a, M* and <j>* in the 
lowest-redshift slice. At higher redshifts we hold a fixed and al- 
low only M* and <j>* to vary. The results of these fits are shown as 
dashed lines in Fig. [T6] and the insets show 95 per cent likelihood 
contours of (M*, lg </>*). These least-squares fits are for illustra- 
tion only. In order to calculate the parameters Qswml and Pswml 



5 This choice of magnitude limits also corresponds closely to those of 
[Blanton et aT||2003b} . 



4.4.2 g-band evolution 

The parametric model provides a good fit out to redshift z ~ 0.2, 
beyond which a bright-end excess results in an over-prediction of 
the LF relative to the non-parametric estimates. From these latter 
estimates, one sees that the number density of blue galaxies is grad- 
ually increasing with redshift, whereas red galaxies show the oppo- 
site trend. Photometric errors are less likely to be a problem in the g 
than in the u band, and so again we suspect that AGNs are dominat- 
ing the bright end of the LF at higher redshifts. This interpretation 
does not necessarily imply rapid evolution of the AGN population 
— the volume sampled at low redshifts is simply too small to detect 
them in significant numbers. 



4.4.3 riz-band evolution 

Evolution in the r, i and z bands is qualitatively very similar, and 
so we discuss them together. The parametric model provides a rea- 
sonable fit in all redshift slices, although it should be said that the 
formal fit probabilities from the likelihood ratio test are mostly 
below 1 per cent. This does not necessarily mean that the model 
is a poor fit, as we see from simulations (Appendix [A} that the 
non-parametric LF estimates are biased when the underlying LF 
is evolving. The likelihood ratio test provides improved probabili- 
ties when we consider narrower ranges in redshift. From the non- 
parametric fits, we see that M* brightens with redshift for all sam- 
ples. At low redshifts, red galaxies have a much higher space den- 
sity than blue, but as redshift increases, the density of red galaxies 
drops and that of blue galaxies increases, until blue galaxies come 
to dominate by redshifts z ~ 0.4. 
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Figure 16. Evolution of the ugriz LFs. The five columns show the ugriz LFs respectively from left to right. The four rows show the LFs in four redshift 
ranges increasing from top to bottom as indicated in the leftmost panels. Filled black squares show SWML estimates for combined red and blue samples, filled 
blue circles and red triangles show SWML LFs for the blue and red samples respectively. Open symbols show the corresponding 1/ Vm ax estimates — in most 
cases these are indistinguishable from the SWML estimates. Continuous lines show the parametric evolving LF for each sample The dotted lines reproduce the 
parametric LF fit for each sample from the lowest redshift bin. Dashed lines show least-squares fits to the SWML estimates with a fixed at higher redshifts. 
The insets show the 95 per cent likelihood contours for (M* , lg <f>* ) parameters obtained from these fits. LFs (but not contours) for the blue and red sample 
have been scaled by a factor of 0. 1 to aid legibility. 
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Table 5. Evolving Schechter function fits to ugriz LFs. Columns are the same as in Table with the addition of. evolution parameters Qpar and P par 
determined from the parametric model and Qswml ar| d J'swml determined from least-squares fits to SWML estimates in eight redshift slices as described 
in the text. 
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Figure 17. Comparison of of our ugriz LFs for the combined sample in three redshift ranges as indicated with previous estimates. Filled and open symbols 
and black lines show our SWML, 1/Vmax and parametric fits respectively. Blue lines (upper panels only) show the z ~ 0.1 SDSS LFs estimated by |Blanton| 
|eFaT]j2003b) . Red lines show the VVDS LFs in bands UBVRI respectively from Ilbert et al. 2005. Green lines show the zCOSMOS S-band LFs from Fig. 1 
of Zucca et al. 2009. The VVDS and zCOSMOS fits are shown as solid lines over the magnitude range actually fitted; dashed lines show extrapolations outside 
the fitted magnitude range. 



4.4.4 Comparison with previous results 

We compare our evolving LFs with previous estimates in Fig. [TTJ 
Here we plot the LFs in three redshift ranges (0.002, 0.2], (0.2, 0.4] 
and (0.4, 0.6]. We choose these ranges to coincide with the first 
three redshift bins used by |Ilbert et al.| p005 ) in their analysis of the 



VIMOS-VLT deep survey (VVDS). Following these authors, we 
assume approximate correspondence between the restframe BVRI 
and the 0A griz passbands, and assume that 0,I) MJ = 0,1 M„ — 
0.25. Their LFs are shown as red lines: solid over the magnitude 
range actually fitted, dashed where extrapolated. 

In the low-redshift range (top row), we also show the non- 
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parametric LF estimates of SDSS galaxies from Blanton et al. 
(2003b) as blue lines. Their estimates are in good agreement with 
our non-parametric SWML and 1/Vmax estimates except that we 
see a slightly higher density of luminous galaxies, particularly in 
the g band ( 0,1 M g — 5 lg ft < — 21 mag). This difference is likely 
to be due to the greater depth of the GAMA sample compared with 
SDSS: the mean redshift in this range is z = 0.13 for GAMA 
versus z = 0.10 for SDSS. The GAMA sample thus contains a 
higher fraction of more distant and hence more evolved galaxies. 
The VVDS Schechter fits in the lowest redshift bin show a slightly 
lower density of luminous galaxies than seen in GAMA and SDSS. 
Note, however, that the bright end of the low-redshift LF is very 
poorly constrained by Ilbert et al. due to their bright apparent mag- 
nitude limit of Iab = 17.5. The zCOSMOS B band LF from Fig. 1 
of |Zucca et al.| ([2009) (Zucca private communication, green line) 
shows a bright-end excess relative to GAMA and to VVDS. While 
reaching about 1.5 mag brighter than VVDS, the zCOSMOS low- 
redshift LF still relies on extrapolation at magnitudes brighter than 
M B - 51gft = -21 mag. 

At intermediate redshifts (middle row), our LFs are in good 
agreement with VVDS and zCOSMOS apart from an excess of it- 
bright galaxies in GAMA, and, conversely, a much higher bright- 
end fit by VVDS in the / ( 01 z) band. This latter discrepancy is 
almost certainly due to poor coverage of the bright end of the / 
band LF at redshifts 0.2 < z < 0.4 by VVDS. 

In the highest redshift range (bottom row), the extrapolation 
of the VVDS Schechter fit shows a higher abundance of luminous 
galaxies in the redder VRI bands than GAMA. In this redshift 
range the comparison is not quite fair, since GAMA contains very 
few galaxies beyond z — 0.5, and so much of the VVDS excess is 
likely to be due to galaxies in the redshift range 0.5 < z < 0.6. 

Overall, our evolving LF estimates are in reasonable agree- 
ment with the previous results o f|Blanton et a l. ( 2003b I, Ilbe rt et al.| 
((2005) and |Zucca et 11^2009) . 



4.4.5 Schechter parameter likelihood contours 

2cr likelihood contours for the Schechter parameters (a, M*), 
determined using the parametric fits, are shown in Fig. [T8] For 
the combined sample, we find roughly consistent faint-end slopes 
a = —1.1 ± 0.1 in all bands. The slightly steeper slope seen in 
the r band is most likely due to the fact that GAMA is selected in 
the r band, and hence we probe slightly further down the LF in this 
than the surrounding bands. Similarly, the slightly steeper slopes in 
our fits to the low-redshift sample (Table [3) are due to the inclu- 
sion of fainter-magnitude galaxies in those fits. As expected, M* 
increases systematically in brightness with increasing wavelength 
from u to z. Red galaxies have systematically shallower faint-end 
slopes (a ~ —0.5 in griz) than blue galaxies (a = —1.45 ± 0.05) 
in all bands. The characteristic magnitudes M* are fainter for red 
galaxies than for blue in the ugr bands and are comparable in the 
iz bands. 

Fig. [18] also shows Schechter function parameters estimated 
from the SDSS main galaxy sample by Blant on et al.| (2003bl. 
We find systematically steeper faint-end slopes (apart from in the 
z band) and brighter characteristic magnitudes (apart from the u 
band). Since our non-parametric estimates are in good agreement 
(Fig. |17) , these differences most likely arise due to strong degen- 
eracies between the parameters a, M* and Q: (a, M*), (a, Q) and 
(M*, Q) are all positively correlated. 
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Figure 18. 2-a likelihood contours for evolving Schechter function param- 
eters a and M* in ugriz bands for combined, blue and red samples (black, 
blue and red contours respectively). Asterisks and dotted ellipses show the 
best-fit values and 2-a error ellipses on the parameters reported by |Blanton] 
|et al.|(2003b) (combined colours only). 
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Figure 19. 2-cr likelihood contours for evolving Schechter function param- 
eters Q and P in ugriz bands for combined, blue and red samples (black, 
blue and red contours respectively). Error bars with symbols show evolu- 
tion parameters and their 2-cr errors determined from a least-squares fit to 
SWML estimates of the LF in eight redshift ranges covering (0.002 < z < 
0.5). Asterisks and dotted ellipses show the best-fit values and 2-cr error 
ellipses on the parameters reported by |Blanton et al.H2003b} (combined 
colours only). 



4.4.6 Evolution parameter likelihood contours 

In Fig.[T9]we show 2a likelihood contours for the luminosity evo- 
lution parameters Q and P from our parametric fits, along with es- 
timates of these quantities and their errors from least-squares fits of 
Schechter functions to the SWML estimates made in eight redshift 
ranges. The differences between the estimates of these parameters 
are frequently larger than the formal errors associated with each 
method. This indicates that our assumption of linear evolution of 
M* and lg <f>* with redshift is only approximate. 

For the combined sample, luminosity evolution is least in the 
r band (Q pa r — 0.7), increasing to Q par — 1.6 in the i and z 
bands. Luminosity evolution is even more pronounced in the g and 
u bands (Q par ~ 2.9 and 6.2, respectively), although, as previously 
noted, the parametric model performs very poorly in these bands, 
and so these values are unreliable at best. Luminosity evolution is 
more pronounced for the red galaxy population than the blue. 

Blue galaxies exhibit positive density evolution, (P par > 
in all bands apart from u and g, Pswml > in all bands apart 
from u), whereas red galaxies show negative density evolution, 
both P par and Pswml are negative in all bands. This observation 
is in good qualitative agreement with an analysis of the zCOSMOS 
survey by |Zucca' et al. (2009), who find that both early- and late- 
spectroscopic-type galaxies brighten in M* by ~ 0.5 mag over the 
redshift range z ~ 0.2 to z ~ 0.9, but that <f>* for early types 
decreases by a factor ~ 1.7 over the same redshift range; for late 
types c6* increases by a factor ~ 1.8. 

Density evolution for the combined sample is positive in the 
r band; in other bands P is either negative or consistent with zero, 
compensating for the stronger luminosity evolution in these bands. 
Thus the contrary density evolution of blue and red galaxies largely 
cancels out in the combined sample. For all bands and samples, the 
evolution parameters (Q, P) are strongly anticorrelated. We remind 
the reader that the maximum-likelihood luminosity evolution Q par 
is determined along with a and M* , independently of the normal- 
isation of the LF. Density evolution P par does depend on the fitted 
value of Q par , as well as the Schechter parameters, resulting in the 
observed anti-correlation between Q and P. In the redder bands, 
riz, the combined LD evolution P + Q is stronger for blue galax- 
ies, (P + Q) P ar ^ 3.7±0.8, than for red, (P + Q) pa r ^ 0.2±0.5. 

Fig. [19] also shows evolution parameters determined from the 
SDSS main galaxy sample by Bl anton et al.| ( |2003b) . The 2a likeli- 
hood contours intersect in gri, and narrowly miss in u and z. In the 
r band we find weaker luminosity evolution and a compensating 
stronger density evolution, vice versa in 2. Our least-squares fits to 
the SWML estimates in the u band yield comparable Q estimates 
to |Blanton et al.| ( |2003b) . Their density evolution, unlike ours, is 
positive, but has a very large error. 

Although sampling a smaller volume, the GAMA data anal- 
ysed here have a mean redshift z ~ 0.2 compared with z ~ 0.1 for 
the data analysed by Blanton et al. We thus have a longer redshift 
baseline over which to measure evolution. 

IHopkins] < [2004> , in an analysis combining constraints from 
the star formation rate density of the Universe and 1.4-GHz ra- 
dio source counts, found Q = 2.70 ± 0.60, P = 0.15 ± 0.60 
for the star forming galaxy population^ This measurement, sensi- 
tive to the star-forming population up to z ~ 1, is consistent with 
our parametric fit results for blue galaxies in the g band at the low- 
redshift end of this range. However, given the very large discrep- 



6 Note that Hopkins actually models evolution as L oc (1 

<$> oc (1 + z) p . 



z)Q and 
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ancy between Qpar and Qswml for blue galaxies in the g band, 
the apparent agreement may be fortuitous. 

For red galaxies, [Brown et al.| (2007) find that M* in the B 
band brightens by ~ —0.7 mag from redshift z = 0.2 to z = 1 
while <f>* declines by about 25 per cent, in qualitative agreement 
with our results 

4.4.7 Evolution summary 

To summarize our findings regarding evolution of the LF: 

(i) The evolutionary model (equation|5j provides a reasonable fit 
in the redder bands, riz, but performs poorly in the u and g bands, 
overpredicting the LF of luminous galaxies at high redshift. This is 
possibly due to a significant contribution from AGNs. 

(ii) Our non-parametric LF estimates are in good agreement 
with SDSS measurements at low redshift and with results from the 
VVDS and zCOSMOS surveys at higher redshifts, over magnitude 
ranges where our LF estimates overlap. 

(iii) There is a strong degeneracy between the luminosity and 
density evolution parameters Q and P. One should be wary in using 
them in isolation, e.g. using the Q parameter to apply evolutionary 
corrections. 

(iv) Nevertheless, red galaxies in all bands show evidence for 
positive luminosity evolution (Q > 0) and negative density evolu- 
tion (P < o). 

(v) Blue galaxies show less luminosity evolution but show evi- 
dence for positive density evolution. 

(vi) The observation of decreasing number density of blue 
galaxies but increasing number density of red galaxies with cosmic 
time implies that the transition from blue cloud to red sequence is 
an important and ongoing phenomenon since redshifts z ~ 0.5. 

(vii) The combined luminosity plus density evolution is stronger 
for blue than for red galaxies. 

4.5 Luminosity density evolution 

As we have seen in the previous section, while it can be difficult to 
isolate the effects of luminosity and density evolution, evolution in 
LD is better constrained. Fig.[20]shows the LD pL sum measured in 
eight redshift bins up to z — 0.5, according to equation[T6] along 
with the prediction pL flt of the parametric model (equation 1 14). 
These results are tabulated in Table|6] For the combined sample, in 
all bands other than it, we see LD increasing with redshift, steeply 
between redshifts z = and z ~ 0.15, slightly more gradually 
thereafter. The blue galaxy LD increases more steeply with redshift 
than the combined sample. The LD of red galaxies barely evolves 
with redshift beyond z ~ 0.15, thus the relative contribution to 
LD from blue galaxies comes to dominate by redshifts z ~ 0.2. 
Given our choice of colour cut (equation]?}, red and blue galaxies 
contribute roughly equally to the LD in the r and i bands at low red- 
shifts (z < 0.15), red galaxies are slightly dominant in the z band 
but under-represented in u and g bands. 

Fig.[20|also s hows the LD estimated from the SDSS by |Blan-| 
Iton et al. (2003b) and |Montero-Dorta & Pradalp009) at a mean 
redshift of z ~ 0.1. Our results are in excellent agreement with 
those of |Blanton et al.| ( |2003b) . |Montero-Dorta & Prada| l |2009} ap- 
pear to have significantly underestimated the <?-band LD, although 
their estimates in other bands are in agreement with ours and those 
of Bla nton et al.| ( |2003b| l. Open triangles show LD estimates from 
CNOC2. We have taken the qo — 0.1 'Total' values from table 3 
of Lin et al. 1999 in U, B and R bands, corresponding roughly 




0.0 0.1 0.2 0.3 0.4 0.5 



z 

Figure 20. Luminosity density in the rive bands ugriz as labelled as a func- 
tion of redshift. Points with error bars show the LD estimated by summing 
galaxies in each luminosity range (equation |16) . Lines show the predictions 
from the parametric fits (equation |14) and shaded regions show confidence 
limits obtained by combining the lower and upper 1-cr limits on Pl(0), 
Q and P. Black squares, blue circles and red triangles show results for 
the combined, blue and red samples respectively. Open circles and squares 
show the LD estimated at redshift z ~ 0.1 from SDSS data by |Blanton] 
|et al.H2003b) and|Montero-Do rta & PradaH2009") respectively. Open trian- 
gles show LD estimates from CNOC2 (Table 3 of Lin et al. 1999) in U, 
B and R bands corresponding roughly to 01 u — 0.25, °' 1 g and 01 i re- 
spectively. Diamonds show u-band LD estimates from Prescott, Baldry & 
|James|(27J09) . 
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Table 6. Luminosity density evolution. The column labelled 'Fit' gives the redshift-zero LD from the parametric fit (equation [l5j, subsequent columns show 
the LD obtained from summing over galaxies (equation] 16) in the indicated redshift ranges. Units are 10 8 Lq/iMpc -3 . 
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01 g and 0-1 i respectively. We convert the CNOC2 



luminosity densities from physical units of h W Hz~ Mpc to 
AB magnitudes using Mab = 34.1 — 2.51gpL and then convert 
into Solar luminosities using the assumed absolute magnitudes of 
the Sun quoted in Section [7!2] The B- and i?-band LD estimates 
are in reasonable agreement with our g-and i-band estimates, while 
the CNOC2 [/-band LD is many times larger than our it-band LD 
at redshifts z > 0.15. Open diamonds show u-band LD estimates 
from |Prescott, Baldry & James| p009) . 

Note that our tt-and <?-band LD estimates will be adversely af- 
fected by the poor fit of the parametric model for LF evolution. The 
selection function in the denominator of equation j 1 6[ i will be over- 
estimated at high redshifts in these bands, and hence the summed 
LD itself will be underestimated, leading to the decline in LD with 
redshift seen for the u band in Fig.|20| 

Previous studies of LD evolution (e.g. |Lilly et al.|1996| 



Lin et| 



al.|1999||Belletal.|20 04; Baldry et al.|2005[|Willmer et al.|2006| 
Faber et al.|2007||Prescott, Baldry & James|2009^ have found that 
for blue and non-colour-selected galaxies, pL increases monotoni- 
cally with redshift, while for red galaxies, it is approximately con- 
stant, with a possible decline beyond redshift z ~ 1. We have pre- 
sented the most detailed investigation to date of LD evolution since 
redshift z ~ 0.5, finding consistent results in the gri and z bands 
with previous analyses which have focused primarily on evolution 
beyond redshifts z ~ 0.5. 



5 CONCLUSIONS 

We have presented the first measurements of the ugriz galaxy 
LFs from the GAMA survey, after correcting for imaging, target 
and spectroscopic incompleteness. At low redshift (z < 0.1), the 
shapes of the blue galaxy LFs are reasonably matched, albeit not in 
detail, by standard Schechter functions. LFs for red galaxies show a 
noticeable dip at intermediate magnitudes, requiring double power- 
law Schechter functions to obtain an adequate fit. One should be 
cautious in interpreting this as the upturn predicted by halo occu- 
pation distribution models (e.g. [Brown et al. 2008} and the |Peng et] 
[aX] ( |2010) quenching model, since the faint end of our red galaxy 
LF contains a significant fraction of edge-on disc systems, which 



are likely to be dust-reddened. We find consistent faint-end slopes 
in all bands, a + f3 = —1.35 ± 0.05 for the combined sample. 

In order to determine evolution of the LF, we employ the para- 
metric model o f|Lin et aT| (T999 ) in which characteristic magnitude 
M* and log density lg (f>* are allowed to vary linearly with red- 
shift. We test the parametric model by comparing with estimates 
using the 1/Knax and SWML estimates. We find that the r, i and 
z bands are qualitatively well fitted by this model, although the 
model provides poor likelihood fits compared with SWML. The 
model predicts an excessively high number density in the u and 
g bands at high redshift, most likely due to QSO/Seyfert contam- 
ination ( Montero -Dorta & P rada 2009). With this caveat in mind, 
we find positive (i.e. increasing with redshift) luminosity evolution 
in all bands and for all colour samples. Luminosity evolution is 
stronger for red than for blue galaxies, with blue galaxies brighten- 
ing by ~ 1-1.5 mag per unit redshift, red galaxies brightening by 
~ 2-2.5 mag per unit redshift. 

Number density evolution for blue galaxies is positive in the 
redder riz bands in which it can be reliably measured, while red 
galaxies exhibit negative density evolution. This observation of de- 
creasing number density of blue galaxies but increasing number 
density of red galaxies with cosmic time implies that the transition 
from blue cloud to red sequence is an important and ongoing phe- 
nomenon since redshifts z ~ 0.5. Investigation of the mechanism 
that causes this transition will be the subject of future work, but it 
appears unlikely that mergers play a dominant role at these moder- 
ate redshifts, given the low merger fraction (~ 5 per cent or less) 
observed at low redshift by e.g. Cons elice et al.| ( [2009| > and |Lotz et] 
[aLl pOTT> . 

Luminosity density increases from redshift zero until z ~ 
0.15, beyond which redshift it increases more gradually for the 
combined sample. The LD of red galaxies is roughly constant be- 
yond z ~ 0.15, whereas that for blue galaxies keeps on increasing, 
leading to blue galaxies dominating the LD at higher redshifts. 

In this paper, we have not considered the effects of internal 
dust extinction on the LF, nor have we considered the effects of 
using total as opposed to Petrosian magnitudes (Graha m et al.|2005| 
Hill et al. 201 1 1. These extensions to the analysis will be considered 
in a future paper, along with a measurement of the galaxy LF for 
AGN-dominated, star-forming and quiescent galaxies which, it is 
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hoped, will resolve the problems encountered while attempting to 
fit an evolutionary model in the u and g bands. 
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APPENDIX A: TESTING THE METHODS 

We here test our methods using simulated mock catalogues. We 
generate clustered distributions of points with known evolving LF 
and then apply the GAMA selection effects to create a set of mock 
catalogues. 



Al Clustered simulations 

We use a Soneira & Peebles] ( |1978[ > type simulation to generate 
a clustered distribution of points in a cube 1200 /i _1 Mpc on a 
side with periodic boundary conditions. The parameters in the sim- 
ulation are chosen to yield similar clustering properties to those 
measured by |Zehavi et al.| ( |2005^ for a flux-limited sample of 
SDSS galaxies, namely a correlation function with power-law slope 
7 ~ 1.7 and correlation length ro — 5.2/i _1 Mpc. This allows us to 
investigate the effects of large-scale structure on our LF estimates. 

With the observer located at one corner of the cube, each 
galaxy in the simulation was assigned an absolute magnitude M 
within the range —24 < M r < — 10 mag, drawn at random 
from a LF with parameters specified in Table |AT| Since these are 
static simulations, we assume a linear redshift-distance relation, 
with r = czj Ho and volume element dV oc r 2 . Apparent magni- 
tudes m are calculated using m = M + 25 + 5 lg r(l + z) with 
no if -correction. Strictly, of course, there is no factor (1 + z) in 
luminosity distance in Euclidean space. We choose, however, to in- 
clude this factor in the simulations to make them more realistic — 
without it one predicts far too many galaxies at higher redshifts. 
The number of observable galaxies -/V b s to redshift z max = 0.4 in 



each simulation cube was determined by integrating the model LF 
4>(L, z) over luminosity and redshift: 

N ohs = / 4>{L,z)dL — dz, (Al) 

JO il mia ( Z ,m lim ) dz 

with mii m = 19.8 mag. 

We crudely mimic the GAMA survey geometry by selecting 
galaxies within each of three 12 x 4 deg 2 regions (bounded by 
latitude 0° < (p < 4° and longitudes 0° < 9 < 12°, 36° < 9 < 
48° and 78° < 9 < 90°). 

Imaging completeness is determined by linear interpolation 
of the curve in Fig[T] r-band half-light surface brightness, /iso,r, 
for each simulated galaxy is assigned according to the empirically 
observed relation between /iso.r and r-band absolute magnitude 
°- 1 M r for GAMA galaxies, fi 50 , r = 22.42 + 0.029° A M r with 
1-sigma scatter of 0.76 mag. Target completeness is obtained using 
the empirically observed completeness shown in Fig. [2] Note that 
we do not attempt to follow the dependence of target complete- 
ness on sky coordinates, and so therefore any dependence of target 
completeness on target density in the real data will not be simu- 
lated. Given the better than 98 per cent completeness of GAMA 
spectroscopy, this should not be a significant issue. Finally, spec- 
troscopic completeness is obtained by generating a fibre magnitude 
for each simulated galaxy according to the empirically observed re- 
lation between Petrosian and fibre magnitudes for GAMA galaxies, 
ffib = 5.84 + 0.747rp ctro with 1-sigma scatter of 0.31 mag. Red- 
shift completeness is then obtained using the sigmoid function fit 
to redshift success shown in Fig. [3] 

Considering each form of completeness in turn, galaxies are 
selected at random with a probability equal to the completeness. 
Weights for the simulated galaxies that survive this culling process 
are assigned in the same way as for observed galaxies (equation|2]l. 

A2 Simulation results 

Eight independent mock catalogues were generated for each of two 
different input LFs, as described above, and used as input cata- 
logues for our LF estimation code. Naturally, when analysing these 
simulations, we assume a consistent cosmology in calculating dis- 
tances, apparent magnitudes and volumes. 

The Schechter function parameters recovered by the paramet- 
ric LF estimator, for both non-evolving and pure luminosity evolu- 
tion simulations, are given in Table [Al] In both cases we recover 
the true LF parameters with minimal bias, ~ 1<t at worst. 

Fig. |Al| shows the LF recovered in eight redshift slices from 
our evolving mock catalogues. In order to amplify any discrepan- 
cies, all estimates have been normalised by the true LF, obtained 
by substituting the input LF parameters into equation i |2 1 fr . We only 
plot binned estimates when there is at least one galaxy in that mag- 
nitude bin in all eight realisations in order to avoid biasing the mean 
high if only realisations with one or more galaxies are included or 
biasing it low if all realisations are used. This eliminates bins fainter 
than M = —13 mag, for which galaxies are only found in a subset 
of the simulations. 

For the low redshift (z < 0.05) slice, all estimates are in good 
agreement with the true LF. Faintwards of M r ~ —16 mag, the 
scatter in 1 /l/ max estimates starts to increase due to density varia- 
tions induced by large-scale structure. The SWML estimator, being 
insensitive to density variations, has smaller error bars at the faint 
end. 

For the higher redshift slices, we see that l/Vmax and SWML 
estimators give essentially identical results. There is a tendency for 
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Table Al. LF parameters estimated from both non-evolving and evolving simulated data. The recovered values show the mean and standard deviation from 
eight mock catalogues. 





a 


M* - 5 lg h 


</>*//i 3 Mpc- 3 


Q 


P 


True 


-1.20 


-20.80 


0.0100 


0.00 


0.00 


Recovered 


-1.19 ±0.01 


-20.78 ± 0.02 


0.0101 ± 0.0007 


0.04 ±0.06 


-0.03 ±0.30 


True 


-1.20 


-20.80 


0.0100 


2.00 


0.00 


Recovered 


-1.20 ±0.01 


-20.79 ±0.02 


0.0101 ± 0.0004 


2.05 ± 0.05 


-0.07 ± 0.13 



0.00<z<0.05, P, =0.213 



0.05<z<0.10, P, =0.234 



0. 0<z<0.15, P, =0.092 



j |ii m i"iiB Mn ij 



0.15<z<0.20, P, =0.136 



fii 



0.20<z<0.25, P, =0.090 



0.25<z<0.30, P, =0.117 



0.30<z<0.35, P, =0.011 



0.35<z<0.40, P, =0.024 



-24 -22 -20 -18 -16 -14 

M-5\oeh 



Figure Al. Simulated evolving LF estimates in eight redshift slices as in- 
dicated. All estimates have been normalised by the 'true' underlying LF of 
the evolving simulations; a perfect estimate would like along the horizon- 
tal dotted line. Open symbols show the mean from eight mock catalogues 
determined using the 1 / V ma x estimator; filled symbols are from SWML es- 
timator. Error bars for each come from the scatter between the eight mocks. 
The continuous lines show the best-fit evolving parametric fit, as given by 
equation |21| Pj, gives the likelihood ratio determination of how well the 
the parametric fit describes the SWML observations in that range. 




0.00 0.05 



0.10 0.15 0.20 0.25 



0.30 0.35 



Figure A2. Luminosity density estimated from simulated data in eight 
redshift slices. Symbols with error bars show the mean and rms from 
eight mock catalogues determined using equation \\(>) . The continuous 
line shows the mean parametric prediction of luminosity density evolution 
(equation |14) and the red dashed line shows the expected evolution, given 
the simulation parameters. 



broader redshift slices the discrepancy worsens. For example, with 
Az = 0.1, the likelihood ratio probabilities are below 0.001 for the 
three higher redshift slices. It is likely that our binned estimates of 
the GAMA LF (Section]?) are biased in a similar fashion. Note that 
|Cole| ( [20"TT| > has recently proposed a method for estimating binned 
LFs whilst simultaneously fitting for luminosity and density evolu- 
tion. 

Fig. |A2| shows the LD estimated from the evolving simula- 
tions. The recovered LD, both in redshift bins, equation ( |16 r , and as 
predicted by the parametric fit, equation ( [14) , is in excellent agree- 
ment with the prediction given the simulation parameters. The de- 
creasing errors at higher redshift indicate that sample variance is 
the largest contributing factor to errors in LD for these simulations. 
This is not the case with the observed LD (Section [4.5fr , where the 
dominant source of error, particularly for the u and g bands, is the 
applicability of the evolution model (equation|5]l. 



both of these binned estimates to slightly underestimate the bright 
end of the LF, a consequence of that fact that both binned esti- 
mates, unlike our parametric fit, make the assumption that the LF 
is independent of redshift. Even dividing the simulation into eight 
redshift slices of width Az — 0.05, there is a systematic change of 
AM* =0.1 across each slice in these Q — 2 simulations. With 
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